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Abstract 

Within the framework of the improved isospin dependent quantum molecular dynamics (Im- 



S3 



On 
O 



X 



IQMD) model, the emission of pion in heavy-ion collisions in the region 1 A GeV as a probe of 
nuclear symmetry energy at supra-saturation densities is investigated systematically, in which the 
pion is considered to be mainly produced by the decay of resonances A(1232) and N*(1440). The 
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■ total pion multiplicities and the n~ /n + yields are calculated for selected Skyrme parameters SkP, 
SLy6, Ska and SIII, and also for the cases of different stiffness of symmetry energy with the pa- 



rameter SLy6. Preliminary results compared with the measured data by the FOPI collaboration 
favor a hard symmetry energy of the potential term proportional to (p/po) 7s with 7 S = 2. 
PACS: 25.75.-q, 13.75.Gx, 25.80.Ls 

Keywords: ImlQMD model; pion emission; Skyrme parameters; symmetry energy 

Heavy-ion collisions induced by radioactive beam at intermediate energies play a significant role 
to extract the information of nuclear equation of state (EoS) of isospin asymmetric nuclear matter 
under extreme conditions. Besides nucleonic observables such as rapidity distribution and flow of 
free nucleons and light clusters (such as deuteron, triton and alpha etc.), also mesons emitted from 
the reaction zone can be probes of the hot and dense nuclear matter. The energy per nucleon in the 
isospin asymmetric nuclear matter is usually expressed as E(p,S) = E(p,5 = 0) + E S y m (p)5 2 + 
0(5 2 ) in terms of baryon density p = p n + p p , relative neutron excess 5 = (p n — p p )/(p n + 
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p p ), energy per nucleon in a symmetric nuclear matter E(p,S = 0) and bulk nuclear symmetry 
energy E S ym = \ d §|r^ |<5=o- In general, two different forms have been predicted by some 
microscopical or phenomenological many-body approaches. One is the symmetry energy increases 
monotonically with density, and the other is the symmetry energy increases initially up to a supra- 
saturation density and then decreases at higher densities. Based on recent analysis of experimental 
data associated with transport models, a symmetry energy of the form -Esym(p) ~ 31.6(p/po) 7 
MeV with 7 = 0.69 — 1.05 was extracted for densities between O.lpo and 1.2p [D 12]. The 
symmetry energy at supra-saturation densities can be investigated by analyzing isospin sensitive 
observables in theoretically, such as the neutron/proton ratio of emitted nucleons, tt~~ /tt + , 
and K°/K + [2]. Recently, a very soft symmetry energy at supra-saturation densities was pointed 
out by fitting the FOPI data [3 J using IBUU04 model [J]. With the establishment of high- 
energy radioactive beam facilities in the world, such as the CSR (IMP in Lanzhou, China), FAIR 
(GSI in Darmstadt, Germany), RIKEN (Japan), SPIRAL2 (GANIL in Caen, France) and FRIB 
(MSU, USA) [2], the high-density behavior of the symmetry energy can be studied more detail 
experimentally in the near future. The emission of pion in heavy-ion collisions in the region 1 A 
GeV is especially sensitive as a probe of symmetry energy at supra-saturation densities. Further 
investigations of the pion emissions in the 1 A GeV region are still necessary by improving transport 
models or developing some new approaches. The ImlQMD model has been successfully applied to 
treat heavy- ion fusion reactions near Coulomb barrier [51 El C] - Recently, Zhang et al analyzed 
the neutron-proton spectral double ratios to extract the symmetry energy per nucleon at sub- 
saturation density with a similar model |8j. To investigate the pion emission, we further include 
the inelastic channels in nucleon-nucleon collisions. 

In the ImlQMD model, the time evolutions of the baryons and pions in the system under the 
self-consistently generated mean-field are governed by Hamilton's equations of motion, which read 
as 



Here we omit the shell correction part in the Hamiltonian H as described in Ref. [6]. The 
Hamiltonian of baryons consists of the relativistic energy, the effective interaction potential and 
the momentum dependent part as follows: 



dH 



dpi' 



dH 



(1) 




(2) 



Here the Pj and rrii represent the momentum and the mass of the baryons. 



The effective interaction potential is composed of the Coulomb interaction and the local inter- 
action 

Uint = Ucoul + Ul oc . (3) 

The Coulomb interaction potential is written as 

Ucoui = \ E e -^erf(r l3 /V4L) (4) 

where the Cj is the charged number including protons and charged resonances. The = [r» — r^l 
is the relative distance of two charged particles. 

The local interaction potential is derived directly from the Skyrme energy-density functional 
and expressed as 

U loc = J V loc {p{v))dv. (5) 
The local potential energy-density functional reads 



2p 1 + 7 Po 2p 2p 
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where the p n , p p and p = p n + p p are the neutron, proton and total densities, respectively, and the 
5 = (p n — p p ) / \p n + p p ) is the isospin asymmetry. The coefficients a, (3, 7, g sur , g l s s ° r , g T are related 
to the Skyrme parameters to, ti,t 2 , £3 and xo, £1, £2, £3 [6]. The parameters of the potential part 
in the symmetry energy term are also derived directly from Skyrme energy-density parameters as 

1 1 

O'sym = — g(2Xo + l)toPo> b sym = — — (2X3 + l)^3Po; 
/g ^2/3 

c S ym = ~— ( ) Po /3 [3^i - h{hx 2 + 4)]. (7) 
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The momentum dependent term in the Hamiltonian is the same of the form in Ref. [9] and 

expressed as 

U mom = 6 - E ^[ln(e( Pj - Pj ) 2 + l)] 2 , (8) 



with 



Pl] ~ (4ttL)3/2 6XP 



(9) 



4L 

which does not distinguish between protons and neutrons. Here the L denotes the square of the 
pocket-wave width, which is dependent on the mass number of the nucleus. The parameters 5 and 
e were determined by fitting the real part of the proton-nucleus optical potential as a function of 
incident energy. 



In Table 1 we list the ImlQMD parameters related to several typical Skyrme forces after 
including the momentum dependent interaction. The parameters a, (3 and 7 are redetermined in 
order to reproduce the binding energy (Eb=-1Q MeV) of symmetric nuclear matter at saturation 
density p an d to satisfy the relation dE J | p=po =0 for a given incompressibility. Combined Eq.(7) 
with the kinetic energy part, the symmetry energy per nucleon in the ImlQMD model is given by 

r, , S 1 ^ 2 \ 2/3 P 7 f PV ( P\ 5/3 , s 

E sym {p) = ~— I "ZTf P ) + dsyrn h b sym — + C sym — . (10) 

6 Am \Z J po \Po J \Po / 

More clearly compared with other transport models, the symmetry energy can be expressed as 

The value 7 S = 1 is used in IQMD model [101 [11]. In Fig. 1 we show a comparison of the energy 
per nucleon in symmetric nuclear matter with and without the momentum dependent potentials 
in the left panel and the nuclear symmetry energy in the right panel for different cases of Skyrme 
forces SkP, Sly6, Ska and SIII from Eq. (10), 7 S =0.5 (soft) and 2 (hard) with C sym =32 MeV in 
Eq. (11), and also compared with the form E sym = 31.6(p/po) M MeV (/i=0.5 and //=2) [TJ. 
Analogously to baryons, the Hamiltonian of pions is represented as 



i=l 



= £ (M+ m l + v f oul > ( 12 ) 



where the p i and m n represent the momentum and the mass of the pions. The Coulomb interaction 
is given by 



N B 

ycorn = ^e^ (13) 

where the N n and A^^ is the total number of pions and baryons including charged resonances. 
Thus, the pion propagation in the whole stage is guided essentially by the Coulomb force. The in- 
medium pion potential in the mean field is not considered in the model. However, the inclusion of 
the pion optical potential based on the perturbation expansion of the A-hole model gives negligible 
influence on the transverse momentum distribution |12j . 

The pion is created by the decay of the resonances A(1232) and N*(1440) which are produced in 
inelastic NN scattering. The cross section of direct pion production is very small in the considered 
energies and not included in the model [13]. The reaction channels are given as follows: 



NN^NA, NN^NN*, NN -<->• AA, 

A <B- A%, N* Nn. (14) 
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The cross sections of each channel to produce resonances are parameterized by fitting the data 
calculated with the one-boson exchange model p~T]. In the 1 A GeV region, there are mostly 
A resonances which disintegrate into a 7r and a nucleon, however, the N* yet gives considerable 
contribution to the high energetic pion yield. The energy and momentum dependent decay width 
is used in the calculation [TB] . 

Pion meson in heavy-ion collisions is mainly produced at supra-saturation densities of com- 
pressed nuclear matter larger than the normal density po- The production of pions is influenced 
by the A (1232) and the Fermi motion of baryons in the vicinity of the threshold energies. The 
7r~/7r + ratio is a sensitive probe to extract the high-density behavior of the symmetry energy 
per energy. Shown in Fig. 2 is a comparison of the measured total pion multiplicity and n~ /n + 
yields by the FOPI collaboration in central 197 Au+ 197 Au collisions [3] and the results calculated 
by IQMD model [TU] as well as by the ImlQMD model for Skyrme parameters SkP, SLy6, Ska 
and SIII, which correspond to different modulus of incompressibility as listed in table 1. The 
total multiplicity of pion is mainly determined by the cross sections of the channels NN -h- iVA. 
The ImlQMD model with four Skyrme parameters predicts rather well the total yields at higher 
incident energies, but slightly overestimates the values near threshold energies, which may be 
influenced by the in-medium cross sections. In this work, we use the in-vacuum cross sections of 
nucleon-nucleon elastic and inelastic collisions. Reasonable consideration of the in-medium inelas- 
tic collisions in producing A and TV* is still an open problem in transport models, which have been 
performed in Giessen-BUU model [TB]. Using the isobar model, one gets the ratio 7r~/7r + =1.95 for 
pions from the A resonance, and 7r~/7r + =1.7 from the N* for the system 197 Au+ 197 Au [T7]. These 
relations are globally valid, i.e. independent of the pion energy. On the other hand, the statistical 
model predicts that the 7t~/tt + ratio is sensitive to the difference in the chemical potentials of 
neutrons and protons by the relation n~ /n + oc exp[2(/x n — fJ, p )/T] = exp[85E sym (p) /T], where the 
T is nuclear temperature [JH]. The observed energy dependence of the n~ /n + ratio is due to the 
re-scattering and absorption process of pions and nucleons in the mean field of the compressed 
nuclear matter. We use the free absorption cross sections in collisions of pions and nucleons by 
fitting the experimental data. The branch ratio of the charged n and 7r° is determined by the 
Clebsch-Gordan coefficients with the decay of the resonances A(1232) and N*(1440). The 7r~ /tt + 
ratio is sensitive to the stiffness of the symmetry energy at the lower incident energies. The Im- 
lQMD model can predict the decrease trend of the ir~ /tt + ratio with incident energy. While the 
ImlQMD model with different Skyrme parameters gives the same excitation functions of the total 
pion multiplicity owing to the same cross sections in the production of pions and resonances for 



each case, the ir~ /n + yields is different resulting from the symmetry energy. 

The compressed nuclear matter with central density about two times of the normal density is 
formed in heavy-ion collisions in the 1 A GeV region. To extract more information of symmetry 
energy in heavy-ion collisions from the pion production, in Fig. 3 we calculated the time evolution 
of average central density from low to high incident energies and the excitation functions of the 
7r~ /tt + ratios with the force SLy6, but different stiffness of the symmetry energy which corresponds 
to hard (7,5=2), linear (7^=1), soft (7 S =0.5) and supersoft (SIII) ) , and also compared with IQMD 
results [TU] as well as the FOPI data |3J. The ImlQMD model gives larger values of 7t~/tt + 
than the ones calculated by IQMD, which mainly results from the cross section of the channel 
Nty — > A and the larger coefficient C sym . We considered the pion absorption process according 
to the Breit- Wigner formula with the cross section given in Ref . [32] . Our calculations show that 
a stiff symmetry energy is close to experimental data. The results does not support a very soft 
symmetry energy at high-density from analyzing the same experimental data reported in Ref. [I] . 
Situation is different in IBUU04 model, each nucleon in the evolution is enforced by the symmetry 
potential associated with isospin and momentum. Inversely, a transport model reported in Ref. 
[T9] also predicted the larger ratios for stiffer symmetry energy from the analysis of the ir~ /tt + and 
K°/K + yields. The influence of the symmetry energy on pion production in heavy- ion collisions 
is also studied from the distribution of transverse momentum of the total charged pions and the 
ratio tt~/tt + for the cases of stiff and soft symmetry energies as shown in Fig. 4. The 7r~ mesons 
are mostly produced from neutron-neutron collisions, and for a stiff symmetry energy, a wider 
high-density zone is formed in the calculation of the ImlQMD model. The larger 7r~ /tt + ratio is 
also clear in the momentum distribution and the larger errors at the higher transverse momentum 
are resulted from the limited simulation events. 

The final 7t~/tt + ratio with different stiffness of the symmetry energy is shown in Fig. 5 
as a function of N/Z of the systems in the reactions 40 Ca+ 40 Ca, 96 Ru+ 96 Ru, 96 Zr+ 96 Zr and 
197 Au+ 197 Au, and also plotted the ratios of N/Z and (N/Z) 2 as a function of N/Z at incident 
energy 0.4A GeV and 1.5 A GeV, respectively. The FOPI data [3] and the results calculated by 
IQMD model [10] are also given for a comparison. Experimental data and calculations show that 
an increase trend of the 7t~/tt + ratio in realistic heavy ion collisions than that predicted by the 
isobar model is found at near threshold energy 0.4A GeV, especially for the larger N/Z systems. 
The ratio decreases with the incident energy and the value is located between the lines of (N/Z) 2 
and N/Z at incident energy 1.5 A GeV. The phenomena can be explained from the fact that the 
symmetry energy enhances the N/Z ratio in the high-density region at lower incident energy. 



The decrease of the tt~/tt + ratio with the incident energy is mainly owing to the production of 
pions from secondary nucleon-nucleon collisions, such as a neutron converts a proton by producing 
7r~. Subsequent collisions of the energetic proton can convert again to neutron by producing tt + . 
One can see that the stiff symmetry energy is also close to the experimental data. Recently, a 
moderately soft symmetry energy with 7 S ~ 0.9 ± 0.3 was extracted from the analysis of neutron- 
proton elliptic flow of the FOPI/LAND data for the reaction 197 Au+ 197 Au using the UrQMD 
model [20]. Further experimental works associated transport models should be performed in more 
details to get reliable information of the high-density trend of the symmetry energy in heavy-ion 
collisions. 

In summary, the pion production in heavy-ion collisions in the region 1 A GeV is investigated 
systematically by using the ImlQMD model. The total multiplicity of produced pion and the 
tt~/tt + ratio in central collisions are calculated for the selected Skyrme parameters SkP, SLy6, 
Ska, SIII which correspond to different modulus of incompressibility of symmetric nuclear matter 
and different cases of the stiffness of symmetry energy, and compared them with the experimental 
data by the FOPI collaborations as well as IQMD results. The 7i^/tt + excitation functions for 
the reaction Au+ 197 Au and the dependence of the tt~/tt + ratio on N/Z of reaction systems 
at energy 0.4A GeV are compared with the force SLy6, but different stiffness of the symmetry 
energy. Calculations show that a stiffer symmetry energy of the potential term with 7 S = 2 is 
close to the experimental data. 

Acknowledgement s 

This work was supported by the National Natural Science Foundation of China under Grant 
Nos. 10805061 and 10775061, the special foundation of the president fund, the west doctoral 
project of Chinese Academy of Sciences, and major state basic research development program 
under Grant No. 2007CB815000. 

References 

[1] L.W. Chen, CM. Ko, B.A. Li, Phys. Rev. Lett. 94 (2005) 032701; L.W. Chen, V. Greco, 
CM. Ko, B.A. Li, Phys. Rev. Lett. 90 (2003) 162701. 

[2] B.A. Li, L.W. Chen, CM. Ko, Phys. Rep. 464 (2008) 113. 

[3] W. Reisdorf, M. Stockmeier, A. Andronic, et al. (FOPI collaboration), Nucl. Phys. A 781 
(2007) 459. 



7 



[4] Z.G. Xiao, B.A. Li, L.W. Chen, et al, Phys. Rev. Lett. 102 (2009) 062502. 

[5] Z.Q. Feng, F.S. Zhang, G.M. Jin, X. Huang, Nucl. Phys. A 750 (2005) 232. 

[6] Z.Q. Feng, G.M. Jin, F.S. Zhang, Nucl. Phys. A 802 (2008) 91; Z.Q. Feng, G.M. Jin, Phys. 
Rev. C 80 (2009) 037601. 

[7] N. Wang, Z.X. Li, X.Z Wu, et al, Phys. Rev. C 69 (2004) 034608. 

[8] Y.X. Zhang, P. Danielewicz, M. Famiano, et al, Phys. Lett. B 664 (2008) 145. 

[9] J. Aichelin, A. Rosenhauer, G. Peilert, et al, Phys. Rev. Lett. 58 (1987) 1926. 

[10] Ch. Hartnack, R.K. Puri, J. Aichelin, et al, Eur. Phys. J. A 1 (1998) 151. 

[11] L.W. Chen, F.S. Zhang, G.M. Jin, Phys. Rev. C 58 (1998) 2283. 

[12] C. Fuchs, L. Sehn, E. Lehmann, et al., Phys. Rev. C 55 (1997) 411. 

[13] B.A. Li, A.T. Sustich, B. Zhang, CM. Ko, Int. J Mod. Phys. E 10 (2001) 1. 

[14] S. Huber, J. Aichelin, Nucl. Phys. A 573 (1994) 587. 

[15] Z.Q. Feng, G.M. Jin, Chin. Phys. Lett. 26 (2009) 062501. 

[16] A.B. Larionov, W. Cassing, S. Leupold, U. Mosel, Nucl. Phys. A 696 (2001) 747; A.B. 
Larionov, U. Mosel, Nucl. Phys. A 728 (2003) 135. 

[17] R. Stock, Phys. Rep. 135 (1986) 259. 

[18] G.F. Bertsch, Nature 283 (1980) 280; A. Bonasera, G.F. Bertsch, Phys. Lett. B 195 (1987) 
521. 

[19] G. Ferini, T. Gaitanos, M. Colonna, et al., Phys. Rev. Lett. 97 (2006) 202301. 

[20] W. Trautmann, M. Chartier, Y. Leifels, el al.. larXiv:0907.28221 



Table 1: ImlQMD parameters and properties of symmetric nuclear matter for Skyrme effective 
interactions after the inclusion of the momentum dependent interaction with parameters 5=1.57 
MeV and e=500 c 2 /GeV 2 



Parameters 


SkM* 


Ska 


SIII 


SVI 


SkP 


RATP 


SLy6 


a (MeV) 


-325.1 


-179.3 


-128.1 


-123.0 


-357.7 


-250.3 


-296.7 


P (MeV) 


238.3 


71.9 


42.2 


51.6 


286.3 


149.6 


199.3 


7 


1.14 


1.35 


2.14 


2.14 


1.15 


1.19 


1.14 


<w(MeV fm 2 ) 


21.8 


26.5 


18.3 


14.1 


19.5 


25.6 


22.9 


C°(MeV f m 2 ) 


-5.5 


-7.9 


-4.9 


-3.0 


-11.3 


0.0 


-2.7 


9r (MeV) 


5.9 


13.9 


6.4 


1.1 


0.0 


11.0 


9.9 


C sym (MeV) 


30.1 


33.0 


28.2 


27.0 


30.9 


29.3 


32.0 


a sym (MeV) 


62.4 


29.8 


38.9 


42.9 


94.0 


79.3 


130.6 


b sym (MeV) 


-38.3 


-5.9 


-18.4 


-22.0 


-63.5 


-58.2 


-123.7 


c sym (MeV) 


-6.4 


-3.0 


-3.8 


-5.5 


-13.0 


-4.1 


12.8 


Poo (fm" 3 ) 


0.16 


0.155 


0.145 


0.144 


0.162 


0.16 


0.16 


m*Jm 


0.639 


0.51 


0.62 


0.73 


0.77 


0.56 


0.57 


i^oo (MeV) 


215 


262 


353 


366 


200 


239 


230 



Q 



2 

p/p 



2 

P/P„ 



Figure 1: The density dependence of the energy per nucleon in symmetric nuclear matter at 
temperature T=0 MeV with and without the momentum dependent potentials (left panel) and 
comparison of the density dependence of the nuclear symmetry energy for different Skyrme forces 
SkP, Sly6, Ska and SIII, and the symmetry energy E sym = 31.6(p/p ) 7 MeV (the two cases 7=0.5 
and 7=2) taken in Refs. [U [2] (right panel). 




Figure 2: Comparison of calculated pion multiplicity and 7r /h + ratios in central 197 Au+ 197 Au 
collisions with different Skyrme parameters, and compared with IQMD results [10] as well as the 
FOPI data [3]. 
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Figure 3: Evolution of average central density at different incident energies (left panel) and the 
excitation functions of the ir~ /tt + ratios at different stiffness of the symmetry energy (hard, linear, 
soft and supersoft), and compared with IQMD results [10] as well as the FOPI data (3] (right 
panel). 




Figure 4: Distributions of transverse momentum of final tt and 7r + and the ratio 7r /n + for the 
cases of stiff and soft symmetry energies in the reaction 197 Au+ 197 Au at incident energy Ei a b = 
0.4A GeV. 
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Figure 5: The 7r /tt + yields as a function of the neutron over proton N/Z of reaction systems for 
head on collisions at incident energy Ei a b = 0.4 A GeV and 1.5 A GeV, respectively. 
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